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Abstract. We have performed the detailed analysis of 174 high-resolution spectra of FGK dwarfs obtained with 
the ELODIE echelle spectrograph at the Observatoire de Haute-Provence. Abundances of Fe, Si and Ni have been 
determined from equivalent widths under LTE approximation, whereas abundances of Mg have been determined 
under NLTE approximation using equivalent widths of 4 lines and profiles of 5 lines. Spatial velocities with an 
accuracy better than 1 kms~^ , as well as orbits, have been computed for all stars. They have been used to define 
2 subsamples kinematically representative of the thin disk and the thick disk in order to highlight their respective 
properties. A transition occurs at [Fe/H] = —0.3. Stars more metal-rich than this value have a flat distribution 
with Zmax < 1 kpc and aw < 20 kms~^ , and a narrow distribution of [a/Fe] . There exist stars in this metallicity 
regime which cannot belong to the thin disk because of their excentric orbits, neither to the thick disk because 
of their low scale height. Several thin disk stars are identified down to [Fe/H] — —0.80. Their Mg enrichment is 
lower than thick disk stars with the same metallicity. We confirm from a larger sample the results of Feltzing 
et al (200'S) and Bensby et al II2003II showing a decrease of [a/Fe] with [Fe/H] in the thick disk interpreted as 
the signature of the SNIa which have progressively enriched the ISM with iron. However our data suggest that 
the star formation in the thick disk stopped when the enrichment was [Fe/H] — —0.30, [Mg/Fe] = -1-0.20, [Si/Fe] 
= -1-0.17. A vertical gradient in [a/Fe] may exist in the thick disk but should be confirmed with a larger sample. 
Finally we have identified 2 new candidates of the HR1614 moving group. 

Key words. Stars: fundamental parameters - Stars; abundances - Stars: kinematics - Stars: atmospheres - Galaxy: 
stellar content 



1. Introduction 

The behaviour of stellar elemental abundances with age 
and kinematics, in various substructures of the Galaxy, 
puts important constraints to the construction of mod- 
els of chemical and dynamical evolution of our galactic 
system. In this paper we put particular attention to the 
transition between the thin disk and the thick disk and 
to the abundances of the a-elements Mg and Si and the 
iron-peak element Ni. According to the current nucleosyn- 
thetic theory, a-elements are being produced as a result 
of a-capture reaction, taking place in the core of massive 
stars during their explosion as SN II (Burbidge et al ll957(l . 
Fe is produced by both massive SN II and less massive SN 
la. It is expected that 2/3 of the solar abundance of iron is 
produced by explosions of white dwarfs in double systems 
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(SN la), 1/3 by SN II (Timmes et al 1995). If the per- 
centage of massive stars in the earlier Galaxy was higher 
than today, one can predict that the a/Fe ratio is going 
to change over time. A well established observational fact 
shows that in old metal-poor stars of the Galaxy, [a/Fe] 
, in particular [Mg/Fe] , is overabundant relative to Sun's 
value (Wallerstein 1961; Gratton & Sneden 1987; Magain 
1989; Nissen et al 1994; Fuhrmann et al 1995 etc). If the 
a-element overabundance is a typical chemical pattern in 
halo stars in comparison with disk stars, there is a ques- 
tion whether there is a distinction in the Mg behaviour in 
other subsystems of the Galaxy, in particular in the thin 
and thick disks. Such a difference would have important 
consequences on the choice of the most probable scenario 
of formation of the thick disk (collapse, accretion etc.) 
and its timescale. Moreover [a/Fe] could serve as an ad- 
ditional criterion for the deconvolution of thick and thin 
disk populations. It is especially important for the stars 
with metallicity in the interval of overlapping [Fe/H] val- 
ues (-0.7 <[Fe/H] < -0.3). 
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Several groups have recently attacked this problem by 
attempting to make consistent analyses of abundances of 
Mg and other elements in large (> 50) samples of stars 
of the thin and thick disks, using kinematics to distin- 
guish the stellar populations. However, their results are 
ambiguous. Fuhrmann (1998) found the thick disk to be 
significantly older than the thin disk and its data exhibit 
an appreciable constant overabundance of Mgl in the thick 
disk relative to the thin disk. He interprets the observed 
discontinuity in the chemical and kinematical behaviour 
of the thin disk and thick disk populations as a hiatus 
in star formation before the earliest stages of the thin 
disk. Mashonkina's group (Mashonkina & Gehren 2000, 
Mashonkina & Gehren 2001, Mashonkina et al 2003) is in 
agreement with Fuhrmann concerning the discontinuity 
between the thin disk and the thick disk. They observe a 
step-like decrease of Eu/Ba ratio at the thick disk to thin 
disk transition and estimate that the thick disk formed 
from well mixed gas during a short timescale of 0.5 
Gyr. Gratton et al (2000) argue in favour of a constant 
and significant overabundance of Mg with respect to Fe in 
the thick disk, a fast formation of the thick disk, a sudden 
decrease in star formation during the transition between 
the thick and thin disk phases and a formation scenario 
mixing dissipational collapse and accretion on similar time 
scales. Contrary to these studies, Chen et al (2000) do 
not observe in their data any discontinuity between the 
thin disk and the thick disk. Their thick disk stars ex- 
hibit a slightly higher a/Fe ratio than the thin disk, but 
not exceeding -1-0.4. Very recently Bensby et al (2003) and 
Feltzing et al H2()()3|l published a convincing study show- 
ing that the thick disk extents at solar metallicities with 
an abundance trend clearly separated from the thin disk 
despite a large overlap in metallicity. A knee and decrease 
in [Mg/Fe] vs [Fe/H] is interpreted as the signature of con- 
tribution from SNIa to the enrichment of the interstellar 
gas out of which the thick disk stars formed. 

However these studies cannot be easily compared be- 
cause they do not rely on the same definition of the thick 
disk. The thin disk and the thick disk are known to greatly 
overlap in kinematics and metallicity, and the selection of 
a sample representative of each population is not obvious. 
The use of a single criterion (Vrot , Z^ax or eccentricity 
of the orbit, [Fe/H] , ...) is usually not sufficient for this 
task. In the last decade, studies of the thick disk kine- 
matics have been numerous and its velocity distribution 
is now well constrained. This allows to make a rigourous 
classification of the stellar populations by computing the 
probability of each star of a sample to belong to the thin 
disk and the thick disk according to its velocity (U,V,W), 
having first evaluated the selection biases which affect lo- 
cal samples. It is the way we have adopted in this paper to 
investigate the behaviour of Mg and Si (as alpha-element) 
and Ni (as iron-peak element) in the thin and thick disks. 

In the past years, Mg abundances in stars have been 
determined in a wide range of metallicities and temper- 
atures, either in LTE approximation or through detailed 
NLTE calculations. In the visible range, there are about 



10 Mgl lines, but most of them are strong lines with equiv- 
alent width (EW) greater than 200 mA, and for such 
lines NLTE effects can be significant. As shown by several 
authors (Zhao & Gehren 2000; Idiart & Thevenin 2000; 
Shymanskaya et al 2000), departures from LTE are ob- 
served in Mgl lines in stars of various spectral types. We 
have therefore used the NLTE approach to determine Mg 
abundances more exactly. The lines that we have used 
for the other elements have small or moderate EW and 
considering our temperature and metallicity interval, we 
have negleted NLTE corrections. As a matter of fact, in 
the Sun the abundance deviations due to NLTE effects are 
generally small and do not exceed 0.1 dex, as established 
in former calculations: 0.07 dex according to Shchukina & 
Bueno (2001), <0.10 dex according to Gehren et al (2001a, 
2001b) for Fcl lines, and -0.01 dex for Sil lines (Wedemeyer 
2001). Gratton et al (1999) and Thevenin & Idiart (1999) 
have investigated NLTE effects for iron lines in dwarfs and 
giants of different metallicity. In both papers, the dominat- 
ing NLTE effect for Fe is the overionization by ultraviolet 
radiation (UV photons). Gratton et al (1999) found that 
NLTE corrections may be neglected in most cases, includ- 
ing the stars on the main-sequence and red giant branch. 
Thevenin & Idiart (1999) derived metallicity corrections 
of about 0.3 dex for stars with [Fe/H] <~3.0, however for 
stars with [Fc/H] >-l this value does not exceed 0.1 dex. 
Abundances of Fe, Si and Ni were therefore determined in 
LTE approximation. 

This paper presents the determination and analysis of 
kinematical parameters and abundances of Mg, Si and Ni 
for 174 dwarf stars in the domain of — 1.0 <[Fe/H] < -1-0.3. 
This study was carried out using a homogeneous spectral 
material, uniform methods of treatment and NLTE calcu- 
lations for 9 Mgl lines. Two subsamples have been defined 
on the basis of kinematics to be representative of the thin 
disk and the thick disk and are used to highlight the chem- 
ical behaviour of the two stellar populations. 

2. Observations and reduction 

All the spectra used in this paper are extracted from 
the most recent version of the library of stellar spectra 
collected with the ELODIE echelle spectrograph at the 
Observatoire de Haute-Provence by Soubiran et al H1998|l 
and Prugniel & Soubiran H2001|l . The performances of 
the instrument mounted on the 193cm telescope, are de- 
scribed in Baranne et al H199fi|l . ELODIE is a very stable 
instrument, build to monitor radial velocity oscillations of 
stars with exoplanets, at a resolving power of 42 000 in 
the wavelength range A A 3850-6800 AA. Spectrum ex- 
traction, wavelength calibration and radial velocity mea- 
surement have been performed at the telescope with the 
on-line data reduction software while straightening of the 
orders, removing of cosmic ray hits, bad pixels and telluric 
lines were performed as described in Katz et al (1998). 

All the spectra of the library have been parametrized in 
terms of (Toff , logg , [Fe/H] ), either from the literature 
or internally with the TGMET code (Katz et al 1998). 
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This allowed us to select a set of 174 FGK dwarfs and 
subgiants (3.7<logg<4.6) for this study, with metalhcities 
in the range —1 <[Fe/H] < +0.3. Several solar spectra 
taken on the Moon and asteroids were also considered as 
references. The selected spectra have a signal to noise ratio 
(S/N) at 5500 A ranging from 100 to 300. 

The continuum level drawing and equivalent width 
measurements were carried out by us using DECH20 
code (Galazutdinov, 1992). Equivalent widths of lines 
were measured by Gaussian function fitting. Their ac- 
curacy was estimated by comparing our measurements 
on solar spectra to those obtained by other authors. 
The mean difference with Edvardsson et al H1993|l is : 
EWEdv - EWour = -2.3 (fj ^ 2.3) mA for 27 lines of 
Fel, Fell, Sil and Nil in common; and with Reddy et al 
iPTH^ : EWRed-EWour = -1-2 (cr = 2.1) mA for 47 lines 
in common. The comparison is shown in Fig. ^ 




20 40 60 80 100 

EW(our) 



Fig. 1. Comparison of solar EW measured in this study 
with those from the literature (Edvardsson et al 119931 - 
filled circles, Reddy et al 120021 - open circles) 



excitation potentials, calibrated from 45 reference stars 
in common with Alonso et al H1996|l . Blackwell & Lynas- 
Gray (fTHMIl and di Benedetto fT5M|- The zero-point of 
the temperature scale was directly adjusted to the Sun, 
based on 11 solar reflection spectra taken with ELODIE, 
leading to the uncertainty in the zero-point of about IK. 
Judging by the small scatter in our final calibrations and 
Tcfi , the selected combinations are only weakly sensitive 
to effects like rotation, micro- and macroturbulence, non- 
LTE and other. The application range of the line-depth 
method is -0.5 <[Fc/H] < +0.5. 

For the most metal-poor stars, Tes was determined 
earlier (Mishenina & Kovthvukh l2001|l . The Hq line- wing 
fitting was used for stars studed in this work. None of 
the stars from Mishenina & Kovthyukh (2001) had their 
temperature measured by line depth ratio because they 
are too metal poor. In order to testify that the temper- 
ature scales adopted in Mishenina & Kovthyukh l)2001|l 
and Kovtyukh et al (2003) are consistent, we show in Fig. 
13 our adopted temperatures versus those estimated by 
Alonso et al fl996), Blackweh & Lynas-Gray l|1998|l and 
di Benedetto 11998 ) . Our temperature scale is slightly hot- 
ter than their by ~ 20 — 30K, as mentioned in Kovtyukh 
et al (2003), but the dispersion is satisfactory (^ 80K, 
Tab. III). We have checked the 2 outhers HD101177 and 
HD165341. We are certain that Alonso et al's tempera- 
ture for HD101177 (5483K) is in error because recently 
Heiter & Luck (12003!) determined a temperature of 6000K 
similar to ours (5932K), a value confirmed by the Hq, 
profile. HD165341 is a variable, active, rotating spectro- 
scopic binary rending its photometric measurements sus- 
picious. On the basis of IR photometry, Alonso et al 
fl996), BlackweU & Lynas-Gray (TU^ and di Benedetto 
fl998,) determined respectively 4978K, 4983K, 4937K. Our 
determination of 5314K agrees well with that of Zboril & 
Byrne lfT!?M| who find Teff =5260K. 



3. Atmospheric parameters 

In order to perform reliable abundance determinations, it 
is crucial to derive accurate stellar parameters, especially 
effective temperatures Toff . V-K or IR flux methods are 
often advocated as the best (Asplund 2003)) but require 
homogeneous infrared observations which are not avail- 
able for all stars. Recently our group has improved the 
line-depth ratio technique pioneered by Gray (1994), lead- 
ing to high precision Toff for most of our program stars 
(Kovtyukh et al 2003). This method, relying on ratios of 
the measured central depths of lines having very different 
functional dependences on Toff , is independent of inter- 
stellar reddening and takes into account the individual 
characteristics of the star's atmosphere. Briefly, a set of 
105 relations was obtained, the mean random error of a 
single calibration being 60-70 K (40-45 K in the most and 
90-95 K in the least accurate cases). The use of ^70-100 
calibrations per spectrum reduced the uncertainty to 5- 
7 K. These 105 relations were obtained from 92 lines, 45 
with low (x <2.77 eV) and 47 with high (x >4.08 eV) 
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Fig. 2. Comparison of our adopted temperatures with 
those estimated by Alonso et al ()199()|l . Blackwell & 
Lynas-Gray (11998(1 and di Benedetto H1998(l . Crosses indi- 
cate stars with [Fe/H] <-0.5. The two outliers HD101177 
and IID165341 are discussed in the text. IID165341 was 
measured by Alonso, Blackwell and Di Benedetto, that is 
why there are 3 points for the same star. 



Finally, Fig. |21 shows a general lack of correlation be- 
tween the Fe and Mg abundances and Toff that testifies 
to a correct choice of the effective temperatures. 
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Table 1. Mean difference and standard deviation between 
our temperatures and those estimated by Alonso et al 
(TWe i, Blackwell & Lynas-Gray lfTM5|l and di Benedetto 
(|1998;1 for 40 stars in common. 



Reference 


N 




a{K) 


Alonso et al II1996II 


43 


31 


84 


Blackwell & Lynas-Gray 


41 


29 


84 


di Benedetto (IB^ 


41 


23 


78 
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Fig. 3. Abundances of Fe and Mg vs Tcff 



Table 2. Number of stars in common, mean difference 
and standard deviation of parameter comparison with 5 
other studies. 



N 


ATeff 


a 


A log g 


a 


A[Fe/H] 


a 


Rer 


37 


-2 


78 


0.03 


0.17 


-0.07 


0.08 


1 


39 


-28 


82 


0.04 


0.22 


-0.04 


0.10 


2 


18 


10 


75 


0.04 


0.15 


-0.01 


0.09 


3 


17 


-99 


56 


-0.06 


0.13 


-0.08 


0.06 


4 


12 


-15 


74 


0.05 


0.16 


0.04 


0.09 
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* 1 - Edvardsson et al (1993); 2 - Gratton et al (1996); 3 
- Fuhrmann (1998); 4 - Chen et al (2000); 5 - Fulbright 
(2000); 



the wings- fitting for the Mgib lines. For all stars the differ- 
ence does not exceed 0.2 dex (such difference was detected 
only for a few program stars). 

Microturbulent velocities Vt were determined by forc- 
ing the abundances determined from individual Fel lines 
to be independent of equivalent width. Starting with 
Vt = l-1 kms~^ , we varied it until abundances computed 
for Fel Hnes (20 mA< EW < 150 mA) and plotted as a 
function of EW-s showed a zero slope. The precision of 
the microturbulent velocity determination is 0.2 kms~^ . 

The abundance of iron relative to solar one [Fe/H] 
=log(Fe/H)* - log(Fe/H)Q was used as the metallicity pa- 
rameter of a star and was obtained from Fel lines. 

The adopted parameters of the target stars are given 
in Table |S1 To check the adopted model atmosphere pa- 
rameters, we compared our values with those derived by 
Edvardsson et al (1993), Gratton et al (1996) Fuhrmann 
dUnHI), Chen et al [ 20QQ and Fulbright ( 2000). The mean 
differences (other - us) and standard deviations are given 
in Table 121 On the whole, the agreement is rougly good. 



The two most commonly used techniques to determine 
surface gravities are the ionization balance of neutral and 
ionized species and the fundamental relation which ex- 
presses the gravity as a function of mass, temperature and 
bolometric absolute magnitude deduced from parallaxe. 
For the latter method, stellar masses can be estimated 
from evolutionary tracks, but metallicities and eventual a- 
enhancement have to be known first. We therefore chooscd 
to derive logg from ionization balance, a method which 
might be affected by NLTE effects. However, a detailed 
study of surface gravities derived by different procedures 
was preformed by AUende Prieto et al H1999fl who con- 
cluded that astrometric and spectroscopic (iron ionization 
balance) gravities were in good agreement in the metallic- 
ity range —1.0 <[Fe/H] < -1-0.3. We compared our adopted 
surface gravities to those determined astrometrically by 
AUende Prieto et al (|1999|l and obtain a mean difference 
and standard deviation of -0.01 and 0.15 respectively for 
39 stars in common. This is consistent with an accuracy 
of 0.1 dex on our estimated spectroscopic gravities. As an 
additional criterion of the reliable choice of log g we used 



4. Abundance determination 

4.1. Fe, Si, Ni abundances 

Fe, Si, Ni abundances were determined from an LTE anal- 
ysis of equivalent widths using the W1DTH9 code and the 
atmosphere models by Kurucz (1993,) . Appropriate models 
for each star were derived by means of standard interpola- 
tion through Tcff and log g . The model metallicities were 
taken with an accuracy of ±0.25 dex around the TGMET 
first guess. Abundances of the investigated elements Fe, 
Si, Ni were carried out from a differential analysis relative 
to the Sun. Solar abundances of Fe, Si and Ni were cal- 
culated with our solar EWs and the oscillator strengths 
loggf from Kovtyukh & Andricvsky (1999), who derived 
loggf for 565 lines of the 27 chemical elements using solar 
spectra. We determined : log A(Fe) = 7.55, log A(Si) = 
7.55, log A(Ni) = 6.25, where log(H) = 12 (according to 
Grevesse & Sauval (1998) log A(Fe) = 7.50, log A(Si) = 
7.55, log A(Mg) = 7.58, log A(Ni) = 6.25). 

Depending on authors, the solar iron abundance varies 
from 7.44 to 7.64 (Asplund et al 2000; Blackwell et al 1995; 
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Table 3. List of Fe, Si and Ni lines, atomic parameters 
and solar EW. Only available at the CDS 



Gehren et al 2000a; Grevesse & Sauval 1998; Shchukina & 
Bueno 2001, etc). This disagreement is caused by several 
factors: the systems of oscillators strengths, the models 
of solar atmosphere (empirical, hydrodynamical, LTE or 
NLTE assumptions), the solar spectrum LTE or NLTE 
synthesis etc. Using laboratory log gf, Blackwell et al 
(1995) have obtained a solar iron abundance of 7.64. From 
calculated and laboratory log gf, Gehren (2001a) found 
an abundance 7.50 and 7.56 from Fe II lines, and 7.47 
and 7.56 from Fe I lines. The detailed hydrodynamical 
models of the Sun result in the solar value of iron abun- 
dance of 7.44 (Fe I lines) and 7.45 (Fe II lines) (Asplund 
et al 120001 2000). In the same time of the meteoretic value 
of the iron abundance is 7.50 (Grevesse & Sauval 1998). 
The same value was found for the solar atmosphere by 
Shchukina & Bueno (2001) using the NLTE approxima- 
tion and hydrodynamical solar model. The solar silicon 
abundance obtained within the framework of the hydrody- 
namical models (Asplund 2000) is sligtly lower (7.51) than 
the meteoric silicon abundance (7.55) (Grevesse & Sauval 
1998). To reduce the influence of the used systems of os- 
cillators strengths, atmosphere models and instrumental 
errors we used a differential method of abundance deter- 
mination. 

The list of Fe, Si and Ni lines, atomic parameters and 
solar EW are given in Table 13 Relative abundances of Fe, 
Si and Ni are given in Table El 

4.2. NLTE calculations for Mg 

The determination of Mg abundance was carried out 
through detailed NLTE calculations using equivalent 
widths of 4 lines (A A 4730, 5711, 6318, 6319 A) and pro- 
files of 5 lines (A A 4571, 4703, 5172, 5183, 5528 A). 

NLTE abundances of Mg were determined with the 
help of a modified version of the MULTI code of Carlsson 
iPH^ described in Korotin et al. fl999a"1999b). We have 
included in the code opacity sources from the ATLAS9 
program (Kurucz I1992|l . This enables a much more accu- 
rate calculation of the continuum opacity and intensity 
distribution in the UV region which is extremely impor- 
tant in the correct determination of the radiative rates of 
h — f transitions. For b—h transitions include only contin- 
uum opacity and for b — f transitions include continuum 
both opacities and line opacities from ATLAS9. A simul- 
taneous solution of the radiative transfer and statistical 
equilibrium equations has been performed in the approx- 
imation of complete frequency redistribution for all lines. 

We employed the model of magnesium atom consist- 
ing of 97 levels: 84 levels of Mgi, 12 levels of Mgii and a 
ground state of Mgili. They were selected from works of 
Martin & Zalubas (fTMn|l and Biemont & Brault ^E^. A 
detailed structure of the multiplets was ignored and each 



Fig. 4. Grotrian diagram for Mgl. Displayed are the ra- 
diative transitions treated explicity in the non-LTE calcu- 
lations 
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LS multiplet was considered as a single term. The fine 
structure was taken into account only for a level 3s3p^P^, 
since it is closely connected to the most important tran- 
sitions in magnesium atom ( A A 2778 - 2782AA , 3829 - 
3838 AA , 5167 A , 5172 A ). 

Within the described system of the magnesium atom 
levels, we considered the radiative transitions between the 
first 59 levels of Mg i and ground level of Mg ii. Transitions 
between the rest levels were not taken into account and 
they were used only in the equations of particle number 
conservation. 

Only transitions with A < 100000 A were selected for 
the analysis. All 424 b — b transitions were included in the 
linearization procedure. 

Grotrian diagram for the Mgl are displayed in Fig. ^ 
Photoionization cross-sections were mainly taken from 
the Opacity Project (Yan et al I1987|l keeping a detailed 
structure of their frequency dependence, including reso- 
nances. For some important b — f transitions, the cross- 
section structure is extremely complicated making it dif- 
ficult to describe it using only simple approximations like 

Oscillator strengths were selected from the extensive 
compilative catalogue by Hirata & Horaguchi H1994fl . 
Some information was obtained through the Opacity 
Project. As we ignored a multiple structure of all the lev- 
els, the oscillator strengths for each averaged transition 
were calculated as f = 
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After the combined solution of radiative transfer and 
statistical equilibrium equations, the averaged levels have 
been split with respect to multiplet structure, then level 
populations were redistributed proportionally to the sta- 
tistical weights of the corresponding sublevels and finally 
the lines of the interest were studied. 

Electron impact ionization was described using 
Seaton's formula fSeaton iT9~6"2|) . Collision rate between 
ground level Mg I and ground level of Mg ii were approx- 
imated by use of the fits from Voronov H1H97(I . For elec- 
tron impact excitation for all allowed b — b transitions we 
used the van Regemorter I)1H()2|I formula. CoUisional rates 
for the forbidden transitions were calculated by using the 
semiempirical formula (Allen . 1973^) , with a coUisional force 
of = 1. 

Inelastic collisions with hydrogen may play a signifi- 
cant role in the atmospheres of cool stars. We took into ac- 
count this effect with the help of Drawin's formulas l|19()8|l 
offered Steenbock & Holweger H1984|l . with a correction 
factor of 1/3. Nevertheless, as it was shown in our test 
calculation, the influence of such collisions is negligible 
for the stars considered in our paper. In particular, the 
variation of the correction factor from to 1 results in a 
change of the equivalent width less than 0.5%. 

For all the transitions we also took into account such 
broadening parameters of lines as radiative damping. 
Stark effect, van der Waals damping and microturbulent 
velocity. For temperatures, considered by us, the influence 
of Stark effect is small. 

The main influence on the profile was exerted by van 
der Waals broadening. The Unsold's l)1955|l formula which 
is widely used to allow for this effect, is known to yield 
somewhat understimated coefficients Cg. To refine these 
coefficients, we conducted a special comparison of the ob- 
served and calculated line profiles in solar spectrum using 
method Shimanskaya, Mashonkina & SakhibuUin (200QJ. 
The Solar Flux Atlas of Kurucz et al. H1984|l and Kurucz's 
model of the solar atmosphere H1992|l were used. To take 
into account the chromospheric growth of temperature, 
this model was completed by a model of the solar chro- 
mosphere from work of Maltby et al. H1986|l . Nevertheless, 
the influence of a chromosphere on the equivalent width 
of considered lines was insignificant (less than 1%). 

The corrections Alog Cg to the classic Unsold (1955) 
coUisional damping constant derived from solar line wing 
fitting are given in Tabled These values are comparable 
for common lines with results obtained in the paper of 
Shimanskaya, Mashonkina & SakhibuUin (2000). Similar 
values were obtained by Barklem, Piskunov & O'Mara 
l|2()() ()1. For example, their correction factors for the lines 
AA 5172, 5183 AA are Alog C6=0.85. OsciUator strengths 
for observed lines were selected from the compilative cat- 
alogue by Hirata & Horaguchi H1994|l . 

The estimates obtained from the NLTE profile analy- 
sis of the Mgl lines in the solar spectrum give the solar 
magnesium content (Mg/H)=7.57±0.02 which agrees well 
with the determination of Grevesse & Sauval (,1998,1 and 



Table 4. Atomic data for the Mg lines 



Lambda x (eV) loggf ACe AGs (Barklem et al 2000) 



4571.10 


0.00 


-5.79 


1.72 




4702 QQ 


4.35 


-0.52 


0.74 


1.33 


4: ( OU.UO 






1 7^ 




5172.68 


2.71 


-0.38 


0.74 


0.85 


5183.60 


2.71 


-0.16 


0.72 


0.85 


5528.40 


4.35 


-0.49 


0.44 


0.95 


5711.09 


4.35 


-1.72 


0.60 




6318.72 


5.11 


-1.97 


1.19 




6319.24 


5.11 


-2.18 


1.19 




6319.49 


5.11 


-2.61 


1.19 





Shimanskaya, Mashonkina & SakhibuUin H2000() . who de- 
rived (Mg/H)=7.58. 

The comparison the NLTE profiles with those observed 
are presented in Fig. NLTE abundances of Mg are given 
in Table El 

The NLTE effects for the used lines appeared to 
be very small (less than 0.05 dex). The similar results 
were obtained in Shimanskaya, Mashonkina & SakhibuUin 




4729 4730 4731 5710 5711 5712 

1 X 




5170 5175 5180 5185 



Fig. 5. Comparison of synthetic NLTE profiles (contin- 
uous line) of the Mgl lines to the observed spectrum of 
HD224930 (dots) 
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Table 5. Effects of uncertainties of model parameters 
and EW measurements on the derived abundances 



HD004614 5965/4.4/1. 1/-0.24 



El 


Teff 


logg 


Vt 


EW 


[Fe/H] 


Total 


Fel 


0.07 


0.02 


0.04 


0.04 


0.00 


0.09 


Fell 


-0.02 


-0.08 


0.02 


0.04 


0.03 


0.10 


Mgl 


0.05 


0.02 


0.02 


0.03 


0.00 


0.06 


Sil 


0.02 


0.01 


0.02 


0.03 


0.01 


0.05 


Nil 


0.07 


0.00 


0.03 


0.04 


0.00 


0.09 




HD022879 


5972/4.5/1. 1/-0.77 




Fel 


0.05 


0.02 


0.02 


0.08 


0.00 


0.10 


Fell 


-0.01 


-0.09 


0.01 


0.07 


0.02 


0.12 


iVigi 


0.07 


0.03 


0.02 


0.04 


0.00 


n nq 


Sil 


0.02 


-0.01 


0.00 


0.05 


0.00 


0.05 


Nil 


0.07 


0.01 


0.01 


0.06 


0.00 


0.09 




HD117176 


5611/4.0/1.0/-0.03 




Fel 


0.07 


0.03 


0.04 


0.04 


0.01 


0.10 


Fell 


-0.04 


-0.10 


0.02 


0.04 


0.04 


0.12 


Mgl 


0.06 


0.03 


0.03 


0.03 


0.00 


0.08 


Sil 


0.01 


-0.01 


0.02 


0.03 


0.02 


0.05 


Nil 


0.06 


0.00 


0.03 


0.04 


0.01 


0.08 



(1) - AT^ff = -100 K; 

(2) - Alogg = +0.2; 

(3) - AVt=+0.2 kms-i ; 

(4) - A[Fe/H] = -0.25 dex; 

(5) - AEW= ± 2mA ; 

(6) - Total error. 



5. Error analysis 

Average standard deviations of abundances obtained from 
72 - 265 lines of Fe I (the number of used lines differs from 
star to star), 11-24 lines of Si I and 12-17 lines of Ni I are 
0.10, 0.08 and 0.08 respectively. Final errors in abundances 
result mainly from errors in the choice of parameters of 
the model atmospheres and in equivalent width measure- 
ments (gaussian fitting, placement of the continuum) in 
the case of Fe, Si, and Ni or fitting a synthetic spectrum 
in the case of Mg. Table [S] lists the errors obtained when 
changing the atmospheric parameters by AToff = -100 K; 
Alogg =+0.2; A Yt=+0.2 kms"^ ; A[Fe/H] =±0.25 and 
by assuming an uncertainty of ±2 mA in the EW. These 
values were adopted taking into account the intrinsic ac- 
curacy of the atmospheric parameter determination, the 
processing of the spectra and the comparison of our pa- 
rameter definition with those of other authors. This test 
has been performed for 3 stars with different characteris- 
tics. 

As seen in Tabled the total uncertainty reaches 0.12 
dex for iron abundance determined from II species and 
0.10 dex in the case of I species. The standard deviation 
obtained by comparing our [Fe/H] determinations to those 
obtained by other authors (Table ^ shows that we are 
consistent with them at a level lower than 0.10 dex. 

For Mg, Si and Ni, the total uncertainty due to pa- 
rameter and EW errors is 0.08 dex, 0.05 dex and 0.09 dex 
respectively. We have also compared our Mg, Si and Ni 



Table 6. Abundance [Mg/H], [Si/Fe], [Ni/Fe] comparison 
with other authors 



N 


Ams 


a 


Asi 


a 


Ajvi 


a 


Ref* 


18 


0.03 


0.09 










1 


10 


0.04 


0.09 










2 


37 






-0.05 


0.05 


0.01 


0.05 


3 


17 






-0.04 


0.05 


0.02 


0.03 


4 


5 






0.01 


0.08 


-0.05 


0.09 


5 



* 1 - Idiard & Thevenin (2000); 2 - Gratton et al (2003); 3 - 
Edvardsson et al (1993); 4 - Chen et al (2000); 5 - Reddy 
et al (2003) 



abundances with those determined by other authors. For 
[Mg/H], we have carried out the comparison with Idiart & 
Thevenin (2000) and Gratton et al (2003) who have also 
used the NLTE approach. Our [Si/Fe] and [Ni/Fe] deter- 
minations have been compared to those of Edvardsson et 
al ijTUMji . Reddy et al PIO^ and Chen et al (2000). The 
results of comparisons are listed in Table |B| The mean 
differences are lower than 0.05 and the standard devia- 
tions lower than 0.1 dex for all elements proving the good 
consistency of our study with previous works relying on 
similar methods. 

6. Stellar kinematics, metallicity, elemental 
abundances 

All the selected stars are bright and nearby enough to have 
parallaxes and proper motions measured by Hipparcos 
fESA [TW7|l . These astrometric quantities have been com- 
bined with radial veloctics measured on the ELODIE spec- 
tra by cross-correlation (with an accuracy better than 100 
ms^^) to compute the 3 components (U,V,W) of the spa- 
tial velocities with respect to the Sun. Combining the mea- 
surement errors on parallaxes, proper motions and radial 
velocities, the resulting errors on velocities are of the order 
of 1 kms""'^ . 

Our first concern was to analyse the content of our 
sample in terms of stellar populations. As we were inter- 
ested in the characterization of abundance patterns in the 
thin disk and the thick disk, wc have selected our sample to 
span the metallicity range —1.0 <[Fe/H] < +0.3 in order 
to define 2 subsamples representative of these two popu- 
lations. A discrimination of thin disk and thick disk stars 
is possible using the fact that the two disks are known to 
be distinct by their spatial distribution and local density, 
and by their velocity, metallicity and age distributions. 
However we have performed this classification with a pure 
kinematical appoach because velocity distributions are the 
best constrained from observations reported in the liter- 
ature. If the thin disk velocity distribution is well known 
thanks to Hipparcos data, authors do not yet agree on 
the properties of the thick disk. A review of the current 
knowledge of the thick disk is given in Norris H1999|l . There 
is still a controversy between the adherents of a flat and 
dense thick disk, with typical scale height of 800 pc and 



8 



Mishenina et al.: [Mg/Fe] vs kinematics in the disk 



Table 7. Kinematical parameters adopted for the thin 
disk and the thick disk from Soubiran et al II2003II 





thin disk 


thick disk 


Viag (kms ^ ) 


-15 


-51 


(7(7 ( kms""*^ ) 


39 


63 


av (kms~^ ) 


20 


39 


aw (kms^^ ) 


20 


39 



I • *• • • •* 
- • •••• - • • 
• • • • • ^ 



local relative density of 6-7% (Reyle & Robin nUUT\i to 
15% (Soubiran et al I2(J(J3|I . and the adherents of a thick 
disk with a higher scale height, typically 1300 pc and a 
lower local relative density, of the order of 2% (Reid & 
Majewski'1993). Velocity dispersions are generally found 
to span typical values between 30 and 50 kms""'^ , with an 
asymmetric drift between -20 and -80 kms^^ . Recently, 
Soubiran et al ( 2003^) determined the kinematical param- 
eters of the thick disk and the thick disk from an unbiased 
sample of clump giants. These values, listed in Tab.|Z| are 
used in the present study. 

The metallicity distributions of the thin and thick disks 
are also matter of debate. Recent thick disk metallicity de- 
terminations quoted mean values from -0.36 fBell liyy^ll to 
-0.70 (Robin et al lTOMI Gilmore et al fTO^ . passing by 
-0.48 (Soubiran et al lMISIl . Morrison et al (1990)) brought 
to the fore low- metallicity stars (—1.6 <[Fc/H] < —1.0) 
with disk- like kinematics. Chiba & Beers ( 2000.) estimate 
that 30% of the stars with -1.6 <[Fc/H] < -1.0 belong to 
the thick disk population. Similarly, Bensby et al II2003|I 
found in their sample stars with super-solar metallicities 
and thick disk kinematics. It remains unclear whether 
these extreme populations are separate from the thick disk 
or their metal-weak and metal-rich tails. Moreover, due to 
the great overlap of the thin disk and thick disk distribu- 
tions, it is very difficult to estimate where the transition 
occurs. Authors involved in the study of metallicity distri- 
butions arc often concerned with the chemical evolution 
of the Galaxy and generally do not attempt to separate 
the thin and thick disks but rather consider the thick disk 
as an integral part of the disk. It was the case for instance 
in the work of Edvardsson et al H1993|l on the age - metal- 
licity relation. Haywood H2001|l performed a revision of 
the solar neighbourhood metallicity distribution and con- 
cluded that it peaks at [Fe/H] =0 and that only 4% of the 
nearby stars have [Fe/H] < —0.50. These considerations 
made us think that the selection of thin disk and thick 
disk stars on a pure kinematical criterion would be more 
reliable. 

Fig.l^lshows the distribution of the sample in the plane 
metallicity - W velocity. It is clear from this plot that a 
transition occurs at [Fe/H] ~ -0.30 : above this value we 
measure a dispersion of the W velocities of 16 kms^^ 
while below this value the dispersion is 38 kms"^ . These 
values are typical of the thin disk and the thick disk re- 
spectively (see Tab |7|l . They are related to the different 
scale heights of the two populations and testify that our 
sample is indeed a mixture of thin and thick disk stars. 



-1 -0.5 

[Fe/H] 

Fig. 6. Distribution of the 174 stars in the plane metal- 
licity - W velocity. 



In order to compute the probability of each star to be- 
long to the thin disk or the thick disk on the basis of its 
spatial velocity, we need to know the kinematical param- 
eters of the two disks and their proportion in our sample. 
The kinematical parameters are listed in Tab. 13 but the 
proportions in our sample are of course different of the real 
proportion in the solar ncibourhood because our selection 
in the range —1 <[Fe/H] < +0.3 is supposed to have in- 
creased the number of thick disk stars comparatively to 
their relative density in the solar neighbourhood (2% to 
15%). Moreover the ELODIE library was build from vari- 
ous observing programs which may have biased its content 
and consequently our sample. 

To estimate the proportion of thin and thick disk stars 
in our sample, we have applied an algorithm of deconvo- 
lution of multivariate gaussian distributions, on (U,V,W) 
velocities. This algorithm, SEM (Celeux & Diebolt ir^SSjl . 
solves iteratively the maximum likelihood equations, with 
a stochastic step, with no a priori information on the 
mixture of the gaussian distributions. The results are un- 
ambiguous, our sample is consistent with a mixture of 
2 gaussian populations with parameters very similar to 
those listed in Tab. [7| but with relative densities of 75% 
and 25% respectively for the thin disk and the thick disk. 
Accordingly, we have computed the probability of each 
star, with a measured velocity (U,V,W), to belong to the 
thin disk (Pri) and to the thick disk {Pr2) : 



Pri 



F, 



1 



- exp — 0.5[-^ 



P2 



F 

4. 



w 



'Wi 



2tt CTij.aVi<JWi "Ui 
F ^ piFi +P2F2 

where pi = 0.75, p2 — 0.25 are the relative densities of the 
2 populations in our sample and Vi, ajj. , ay. , ,i — 1,2 
are the kinematical parameters of the thin and thick disks 
listed in Table □ 

From these probabilities we have selected two subsam- 
ples kinematically representative of each population in or- 
der to highlight their respective properties. Putting the 
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limit at Pri > 80% ensures a minimal contamination of 
each subsample by the other population. The thin disk 
subsample includes 109 stars whereas the thick disk sub- 
sample include 30 stars. In the figures, unclassified stars 
are represented by small dots, thick disk stars by larger 
black dots, thin disk stars by open dots. The only halo 
star of the sample, HD194598 ([Fe/H] =-1.21, V=-276.40 
kms^^ ) is not represented in the plots for a better clarity. 

Fig. [71 IHl and |^ show the trend of the abundances of 
Mg, Si and Ni as a function of [Fe/H] . It can be seen that 



— there are a few thin disk stars with [Fe/H] < —0.30 

— metal poor stars ([Fe/H] < —0.60) are all enriched in 
Mg and Si ([Mg/Fe] > -1-0.20, [Si/Fe] > +0.15) 

— on the contrary at solar metallicity, the enrichment of 
Mg and Si does not exceed -t-0.20 

— at a given metallicity thick disk stars have higher 
[Mg/Fe] on average than thin disk stars 

— the dispersion of [Si/Fe] is remarkably small for [Fe/H] 
> —0.30 but quite high at lower metallicity 

— [Mg/Fe] and [Si/Fc] decline with metallicity from 
about -hO.40,-1-0.35 to 0.0, -1-0.08 respectively 

— there are stars with thick disk kinematics at solar 
metallicities, their abundance trends follow the thin 
disk 

— a rise and upturn is visible for [Ni/Fe] for [Fe/H] > 

All these features are in very good agreement with 
those described in Bensby et al l|2()()3(l . 
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Fig. 8. Same as Fig.Qfor Si. 
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Fig. 9. Same as Fig. [3for Ni. 



0.5 




[Fe/H] 

Fig. 7. [Mg/Fe] vs [Fe/H] for the whole sample. Black 
dots indicate thick disk stars, open dots thin disk stars, 
small dots represent the unclassified stars. 



Fig. El represents the correlation between Mg and Si 
relative to iron. Mg and Si are a elements which are sup- 
posed to be mainly produced in SNII. It can be seen in this 
plot that a transition occurs at [Mg/Fe] ~ -1-0.2. Stars with 
[Mg/Fe] < +0.2 have a mean abundance of [Si/Fe] =-f 0.08 
with a very low dispersion of ±0.03, lower than our error 
estimates. [Mg/Fc] is more dispersed (±0.06) around the 
mean of [Mg/Fe] =+0.05. For stars with [Mg/Fe] > +0.2 
the distribution is consistent with a linear correlation : 
[Si/Fe] =0.7[Mg/Fe] +0.06 (rms=0.06). 
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Fig. 10. Correlation between [Mg/Fe] and [Si/Fe] 
7. Discussion 

7.1. Stars with thick disk kinematics and thin disk 
metallicity 

A remarkable feature in the Fig.jTlto^lis the behavior of 
stars with thick disk kinematics and thin disk metallicity. 
They follow exactly the chemical trend of the thin disk. 
This led us to wonder if these stars were real thick disk 
stars. According to the W vs [Fe/H] distribution shown 
in Fig. 1^1 metal-rich stars are expected to have a smaller 
scale height than metal poor stars. To have confirmation, 
we have computed the orbital parameters of our stars, by 
integrating the equations of motion in the galactic model 
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of Allen & Santillan H1991|l over an age of 8 Gyr. The 
adopted velocity of the Sun with respect to the LSR is (9, 
5, 6) kms~^ , the solar galactocentric distance Rq = 8.5 
kpc and circular velocity Yisr = 220 kms~^ . Table |H1 
lists the spatial velocities and the parameters of the orbits 
: apogalactic and perigalactic distances (Rmin, Rmax), the 
maximum distance from the plane (Zmax ) and eccentricity 
(e). The distribution of Z^ax vs [Fe/H] in Fig.llllconfirms 
the small scale height of stars with thick disk kinematics 
and solar metallicity, ruling them out of the thick disk. 
They have been classified into the thick disk because of 
they have a large eccentricity, but they have a scale height 
which is typical of the thin disk. Thus the origin of these 
stars cannot be in the thick disk neither in the local thin 
disk despite the same chemical behaviour. The origin of 
these stars has to be found elsewhere. The only ambigu- 
ous star, HD190360 having [Fe/H] =+0.12 and Z^ax =1 
kpc, appears as an outlier in the distribution. This star is 
discussed farther. 
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Fig. 11. The maximum distance to the plane of orbits, 
Zmax , versus [Fe/H] 



to demonstrate if they belong or not to the thick disk we 
have calculated the orbits of their stars and classified them 
with the same criterion as ours (Pri > 80%). The distribu- 
tion of Zniax versus [Fe/H] is shown on Fig. 12. Globally 
their distribution is similar to ours despite their sample 
is much smaller. There is a sudden rise in Z^ax at [Fe/H] 
< —0.3. They have also 3 outliers with solar metallici- 
ties and Z^ax ^ 1 kpc, HD190360 being in common with 
us. The five other stars with thick disk kinematics and 
[Fe/H] > —0.30 have a low scale height typical of the thin 
disk in agreement with our findings. The question wether 
HD190360, HD003648 and HD145148 are really part of 
the thick disk is an open question because they appear as 
outliers in the distribution of kinematics vs metallicity. It 
is only with complete samples with no kinematical bias 
that the relation of such stars with the thick disk will be 
clarified. 
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Fig. 12. The maximum distance to the plane of or- 
bits, Zmax , versus [Fe/H] for Bensby et al's sample. The 
classification thin disk / thick disk is the same as ours 
(Pr, > 80%). 



Among the 8 stars that verify Pr2 > 80% and [Fe/H] 
> -0.30, HD030562, HD139323, HD139341 are K dwarfs 
with very similar velocities and metallicities (Tab. [SJ 
+0.18 <[Fe/H] < +0.21. HD030562 has been identified 
to be part of the HR 1614 moving group by Feltzing & 
Holmberg ( 2000). It is very likely that we have discovered 
2 new candidates of this moving group whose age is esti- 
mated to be 2 Gyr. The origin of this group is however 
unknown. 

HD010145, HD135204, HD152391 have \U\ > 80 
kms^^ . An hypothesis is that they have been thrown 
out from the inner disk by the galactic bar. The over den- 
sity of metal-rich stars with excentric orbits confined near 
the plane have have been interpreted as the signature of 
the galactic bar in the solar neighbourhood by Raboud et 
al 1,1998 ) for instance. 

Bensby et al (|2()03|l have also found several stars with 
thick disk kinematics and thin disk metallicity on the ba- 
sis of a similar probability calculation than ours, but using 
the ratio of Pr2 to Pri. If in our case these stars have a 
clearly a flatter distribution than the thick disk, it was 
not mentioned in Bensby et al's study. As it is important 



7.2. The thin disk 

Several thin disk stars are found at low metallicity, down 
to [Fe/H] = —0.80. Several previous studies have also es- 
tablished the existence of such stars, however the exact 
contribution of metal-poor stars to the thin disk is difficult 
to evaluate from uncomplete samples. Rcddy ct al f2 002|l 
have found in their samples a significant number of stars 
with [Fc/H] < —0.40 that they identify as belonging to 
the thin disk. They have analysed 181 FG dwarfs, a sam- 
ple size similar to ours. Their sample spans the metallic- 
ity interval -0.70 to +0.15, with a peak at [Fe/H] =-0.40. 
Only 3 stars were identified as belonging to the thick disk, 
based on the criterion [JVIg/Fe] > +0.20, all the other stars 
being identified as thin disk stars. The peak at [Fe/H] =- 
0.40 shows clearly that they have favoured in some way 
the metal-poor tail of the thin disk. Unfortunatly their 
selection criteria are not described. On the contrary, the 
way Chen et al f2000| have favoured metal-poor thin disk 
stars in their sample is very clear : they have selected 
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equal number of stars in every metallicity interval of 0.1 
dex in the range —1.0 <[Fe/H] < +0.3. According to 
the low density of the thick disk in the solar neighbour- 
hood, thin disk stars are still dominating at moderatly 
low metallicity. Haywood H2001fl has very carefully taken 
into account the selection biases of his sample to study the 
metallicity distribution of the local disk but his metallic- 
ities are photometric and no special attention was paid 
to a kinematical distinction between the thin disk and 
the thick disk. Nevertheless he found, contrary to Wyse 
& Gilmorc (1995|), a low contribution of the thin disk at 
[Fe/H] < -0.40. 

Our observations suggest that the distribution of 
[a/Fe] in the thin disk is very narrow, specially for Si, 
at [Fe/H] > —0.30. On this point we are in perfect agree- 
ment with Reddy et al PUU^ and Bensby et al 1^011^ . For 
the metal-rich part, we obtain as mean values and disper- 
sions : [Mg/Fe] =-f0.05, cr[Mg/i^e]=0.07, [Si/Fe] =+0.07, 
o'[Sj/Fe]=0.03. Such a narrow chemical distribution implies 
that the stars formed from an homogeneous gas. At lower 
metallicity, the Mg enhancement is higher. However at a 
given metallicity the thin disk shows a lower Mg enhance- 
ment than the thick disk. The behaviour of Si is quite 
surprising, showing high dispersion, whereas Bensby et al 
H2()()3|l find a similar behaviour than Mg. We notice also 
that the Sun would be slightly Si-poor as compared to 
the thin disk. However this offset stays within the error 
bars and might be related to the adopted zero point of our 
abundance scales. Many previous studies (Edvardsson et 
allB?^ Chen et al HUM Reddy et ali2002, Bensby et al 
I2()()3|l exhibit also a positive offset in [Si/Fe] of the thin 
disk. 

A rise of [Ni/Fe] from 0.1 at [Fc/H] =-0.20 to - 

-1-0.1 at [Fe/H] =-f-0.20 is visible in Fig. El as well as an 
upturn for the most metal rich star. A similar feature is 
also observed in Bensby et al l|2()()3|l . 



7.3. The thick disk 

Having eliminated from the thick disk the 8 stars having 
high metallicity and a flat distribution fSect. m|l . we see 
from Figure [7| that [Mg/Fe] decreases from the halo value 
i+O.AO) at [Fe/H] =-1.00 to +0.20 at [Fe/H] =-0.30. Si has 
the same behaviour, decreasing from +0.35 to +0.17, with 
one deviating star HD110897 having a lower enrichment 
in Si. These findings are nicely consistent with the trends 
found by Bensby et al H2()()3|l that they interpreted as the 
signature of the chemical enrichment by SNIa. Feltzing et 
al H2003|l claimed the existence of a "knee" in the [a/Fe] 
abundance trends tracing the beginning of the contribu- 
tion of the SNIa. This " knee" may also exist in our data, 
especially in [Mg/Fe] but the dispersion is too high to 
conclude. Definitively a larger sample of thick disk stars 
is necessary to clarify this point. Contrary to Bensby et al 
H2003|l and Feltzing et al H2003|l . our data suggest that the 
metallicity distribution of the thick disk docs not continue 
at solar metallicity. The trend of declining ehancement. 



which was also mentioned in Prochaska et al H200UI has 
to be compared to the results of Fuhrmann H1998|l and 
Gratton et al H2()()()|l who found in the thick disk a con- 
stant overabundance of a-elements with respect to iron, 
with a value similar to that of the halo. In Fuhrmann 
(1998), only thin disk stars exhibit a decrease of [Mg/Fe] 
with increasing [Fe/H] . We believe that we have been able 
to detect this trend, as Feltzing et al H2003|> and Bensby 
et al l|2003|l . because we have used a more precise crite- 
rion to isolate a sample a pure thick disk stars. We agree 
with most of the previous works that a discontinuity ex- 
ists between the chemical distributions of the thin and 
thick disks. Our data suggest that the star formation in 
the thick disk stopped when the gas had the composition 
[Fe/H] =-0.3, [Mg/Fe] =+0.20, [Si/Fe] =+0.17. 

The distribution of Zmax vs [Fe/H] fFig. [TT|l does not 
exhibit a clear vertical gradient in the thick disk metal- 
licity. We have searched for a vertical gradient in the a- 
element abundances. Fig. El represents Z^ax vs [a/Fe] 
=0.5([Mg/Fe] +[Si/Fe] ) for the thick disk (Pr2i0.80 and 
[Fe/H]i-0.3). A gradient may exist, but should be con- 
firmed with a larger sample. The existence of a vertical 
gradient would imply a significant timescale for the for- 
mation of the thick disk. This implication is also true if 
the decline of a enhancement is interpreted by the onset 
of SNIa. 



. HDl 10897 
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Fig. 13. A possible vertical gradient in [a/Fe] in the thick 
disk ? 



8. Summary and perspectives 

We have investigated the chemical behaviour of the thin 
disk and the thick disk with a large sample of FGK dwarfs 
spanning the metallicity interval [-1.00 ; +0.30]. We have 
used an homogeneous set of high resolution, high S/N 
echelle spectra to derive abundances of Fe, Mg, Si and Ni 
with detailed NLTE calculations for Mg. All sources of er- 
rors have been analysed and we have compared our abun- 
dances with several studies taken in the literature, show- 
ing a good agreement. Orbits and space velocities were 
computed and used to estimate the probability of each 
star to belong to the thin disk and the thick disk, taking 
into account the selection bias that affects our local sam- 
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pie where the thick disk is overrepresented as compared 
to the solar neighbourhood. Our main findings concern 
the transition in kinematics and abundance trends which 
occurs at [Fe/H] = —0.30, the decline of a enhancement 
with increasing [Fe/H] in the thick disk, and the narrow 
distribution of Mg and Si in the thin disk. We have also 
observed several stars which cannot belong to the thin disk 
because of their large eccentricity, neither to the thick disk 
because of their low scale height. This population, perhaps 
coming from the inner disk due to a dynamical effect, has 
to be taken into account when studying the metallicity 
distributions of the stellar populations in the solar neigh- 
bourhood. Among these stars we have discovered two new 
candidates of the HR 1614 moving group. 

Larger samples are necessary to confirm our findings 
on the thin disk and the thick disk. Unfortunatly thick 
disk stars are quite rare in the solar neighbourhood and 
the construction of a significant sample of local thick disk 
stars implies either a huge complete sample of disk stars, 
among which 2% to 15% of thick disk stars are expected, 
or the selection of such stars by kinematical or chemical 
criteria which bias the conclusions which can be drawn. 
A lot of observing material is available. In the public 
ELODIE library, there are hundreds of high quality spec- 
tra of FGK stars that have not been analysed yet. Most of 
these stars have accurate parallaxes and proper motions 
from Hipparcos and their radial velocities are also avail- 
able. It is thus possible to compute their probability to 
belong to a given kinematical population. More difficult 
but necessary is to evaluate the biases which render the 
metal-poor and high velocity stars more represented in the 
library than in the solar neighbourhood, reflecting the in- 
terest of observers in this kind of stars. Another solution 
to avoid these biases is to observe the thick disk in situ 
that is at distances and galactic latitudes where it begins 
to dominate the thin disk in density. We are continuing 
such an observing program (Soubiran et al 12003(1 . 

Detailed analysis of spectroscopic data is a tedious 
work that takes much time to measure equivalent widths, 
to fit profiles, to compute models etc... It is now abso- 
lutly necessary to develop automatic methods to take ad- 
vantage of the huge amount of high quality data which 
is provided by the new generation of spectrographs. We 
have undertaken such work to continue our investigation 
of the correlation between elemental abundances and kine- 
matics among galactic disk stars from larger, more distant 
and complete samples. 
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Table 8. Atmospheric parameters, abundances, 3D velocities with respect to the Sun, orbital parameters, probability to belong to the thin disk and the thick 
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